Trypanosoma cruzi mitochondrial maxicircles display species- and strain-specific variation and a conserved element in the non-coding region

Abstract Background The mitochondrial DNA of kinetoplastid flagellates is distinctive in the eukaryotic world due to its massive size, complex form and large sequence content. Comprised of catenated maxicircles that contain rRNA and protein-coding genes and thousands of heterogeneous minicircles encoding small guide RNAs, the kinetoplast network has evolved along with an extreme form of mRNA processing in the form of uridine insertion and deletion RNA editing. Many maxicircle-encoded mRNAs cannot be translated without this post-transcriptional sequence modification. Results We present the complete sequence and annotation of the Trypanosoma cruzi maxicircles for the CL Brener and Esmeraldo strains. Gene order is syntenic with Trypanosoma brucei and Leishmania tarentolae maxicircles. The non-coding components have strain-specific repetitive regions and a variable region that is unique for each strain with the exception of a conserved sequence element that may serve as an origin of replication, but shows no sequence identity with L. tarentolae or T. brucei. Alternative assemblies of the variable region demonstrate intra-strain heterogeneity of the maxicircle population. The extent of mRNA editing required for particular genes approximates that seen in T. brucei. Extensively edited genes were more divergent among the genera than non-edited and rRNA genes. Esmeraldo contains a unique 236-bp deletion that removes the 5'-ends of ND4 and CR4 and the intergenic region. Esmeraldo shows additional insertions and deletions outside of areas edited in other species in ND5, MURF1, and MURF2, while CL Brener has a distinct insertion in MURF2. Conclusion The CL Brener and Esmeraldo maxicircles represent two of three previously defined maxicircle clades and promise utility as taxonomic markers. Restoration of the disrupted reading frames might be accomplished by strain-specific RNA editing. Elements in the non-coding region may be important for replication, transcription, and anchoring of the maxicircle within the kinetoplast network.

Unique Alterations of an Ultraconserved Non-Coding Element in the 3 ' UTR of ZIC2 in Holoprosencephaly

Coding region alterations of ZIC2 are the second most common type of mutation in holoprosencephaly (HPE). Here we use several complementary bioinformatic approaches to identify ultraconserved cis-regulatory sequences potentially driving the expression of human ZIC2. We demonstrate that an 804 bp element in the 3' untranslated region (3'UTR) is highly conserved across the evolutionary history of vertebrates from fish to humans. Furthermore, we show that while genetic variation of this element is unexpectedly common among holoprosencephaly subjects (6/528 or >1%), it is not present in control individuals. Two of six proband-unique variants are de novo, supporting their pathogenic involvement in HPE outcomes. These findings support a general recommendation that the identification and analysis of key ultraconserved elements should be incorporated into the genetic risk assessment of holoprosencephaly cases.

Expression of Mitochondrial Non-coding RNAs (ncRNAs) Is Modulated by High Risk Human Papillomavirus (HPV) Oncogenes

The study of RNA and DNA oncogenic viruses has proved invaluable in the discovery of key cellular pathways that are rendered dysfunctional during cancer progression. An example is high risk human papillomavirus (HPV), the etiological agent of cervical cancer. The role of HPV oncogenes in cellular immortalization and transformation has been extensively investigated. We reported the differential expression of a family of human mitochondrial non-coding RNAs (ncRNAs) between normal and cancer cells. Normal cells express a sense mitochondrial ncRNA (SncmtRNA) that seems to be required for cell proliferation and two antisense transcripts (ASncmtRNAs). In contrast, the ASncmtRNAs are down-regulated in cancer cells. To shed some light on the mechanisms that trigger down-regulation of the ASncmtRNAs, we studied human keratinocytes (HFK) immortalized with HPV. Here we show that immortalization of HFK with HPV-16 or 18 causes down-regulation of the ASncmtRNAs and induces the expression of a new sense transcript named SncmtRNA-2. Transduction of HFK with both E6 and E7 is sufficient to induce expression of SncmtRNA-2. Moreover, E2 oncogene is involved in down-regulation of the ASncmtRNAs. Knockdown of E2 in immortalized cells reestablishes in a reversible manner the expression of the ASncmtRNAs, suggesting that endogenous cellular factors(s) could play functions analogous to E2 during non-HPV-induced oncogenesis.

Androgen responsive intronic non-coding RNAs

Abstract Background Transcription of large numbers of non-coding RNAs originating from intronic regions of human genes has been recently reported, but mechanisms governing their biosynthesis and biological functions are largely unknown. In this work, we evaluated the existence of a common mechanism of transcription regulation shared by protein-coding mRNAs and intronic RNAs by measuring the effect of androgen on the transcriptional profile of a prostate cancer cell line. Results Using a custom-built cDNA microarray enriched in intronic transcribed sequences, we found 39 intronic non-coding RNAs for which levels were significantly regulated by androgen exposure. Orientation-specific reverse transcription-PCR indicated that 10 of the 13 were transcribed in the antisense direction. These transcripts are long (0.5–5 kb), unspliced and apparently do not code for proteins. Interestingly, we found that the relative levels of androgen-regulated intronic transcripts could be correlated with the levels of the corresponding protein-coding gene (asGAS6 and asDNAJC3) or with the alternative usage of exons (asKDELR2 and asITGA6) in the corresponding protein-coding transcripts. Binding of the androgen receptor to a putative regulatory region upstream from asMYO5A, an androgen-regulated antisense intronic transcript, was confirmed by chromatin immunoprecipitation. Conclusion Altogether, these results indicate that at least a fraction of naturally transcribed intronic non-coding RNAs may be regulated by common physiological signals such as hormones, and further corroborate the notion that the intronic complement of the transcriptome play functional roles in the human gene-expression program.

Identification of protein-coding and non-coding RNA expression profiles in CD34+ and in stromal cells in refractory anemia with ringed sideroblasts

Abstract Background Myelodysplastic syndromes (MDS) are a group of clonal hematological disorders characterized by ineffective hematopoiesis with morphological evidence of marrow cell dysplasia resulting in peripheral blood cytopenia. Microarray technology has permitted a refined high-throughput mapping of the transcriptional activity in the human genome. Non-coding RNAs (ncRNAs) transcribed from intronic regions of genes are involved in a number of processes related to post-transcriptional control of gene expression, and in the regulation of exon-skipping and intron retention. Characterization of ncRNAs in progenitor cells and stromal cells of MDS patients could be strategic for understanding gene expression regulation in this disease. Methods In this study, gene expression profiles of CD34+ cells of 4 patients with MDS of refractory anemia with ringed sideroblasts (RARS) subgroup and stromal cells of 3 patients with MDS-RARS were compared with healthy individuals using 44 k combined intron-exon oligoarrays, which included probes for exons of protein-coding genes, and for non-coding RNAs transcribed from intronic regions in either the sense or antisense strands. Real-time RT-PCR was performed to confirm the expression levels of selected transcripts. Results In CD34+ cells of MDS-RARS patients, 216 genes were significantly differentially expressed (q-value ≤ 0.01) in comparison to healthy individuals, of which 65 (30%) were non-coding transcripts. In stromal cells of MDS-RARS, 12 genes were significantly differentially expressed (q-value ≤ 0.05) in comparison to healthy individuals, of which 3 (25%) were non-coding transcripts. Conclusions These results demonstrated, for the first time, the differential ncRNA expression profile between MDS-RARS and healthy individuals, in CD34+ cells and stromal cells, suggesting that ncRNAs may play an important role during the development of myelodysplastic syndromes.

Association study between the rs165599 catechol-O-methyltransferase genetic polymorphism and schizophrenia in a Brazilian sample

Schizophrenia is a severe psychiatric disorder with frequent recurrent psychotic relapses and progressive functional impairment. It results from a poorly understood gene-environment interaction. The gene encoding catechol-O-methyltransferase (COMT) is a likely candidate for schizophrenia. Its rs165599 (A/G) polymorphism has been shown to be associated with alteration of COMT gene expression. Therefore, the present study aimed to investigate a possible association between schizophrenia and this polymorphism. The distribution of the alleles and genotypes of this polymorphism was investigated in a Brazilian sample of 245 patients and 834 controls. The genotypic frequencies were in Hardy-Weinberg equilibrium and no statistically significant differences were found between cases and controls when analyzed according to gender or schizophrenia subtypes. There was also no difference in homozygosis between cases and controls. Thus, in the sample studied, there was no evidence of any association between schizophrenia and rs165599 (A/G) polymorphism in the non-coding region 3' of the COMT gene.

Analysis of Hepatitis C Virus Intrahost Diversity across the Coding Region by Ultradeep Pyrosequencing

Hepatitis C virus (HCV) is the leading cause of liver disease worldwide. In this study, we analyzed four treatment-naive patients infected with subtype 1a and performed Roche/454 pyrosequencing across the coding region. We report the presence of low-level drug resistance mutations that would most likely have been missed using conventional sequencing methods. The approach described here is broadly applicable to studies of viral diversity and could help to improve the efficacy of direct-acting antiviral agents (DAA) in the treatment of HCV-infected patients.

Structural studies of Helicase NS3 variants from Hepatitis C virus genotype 3 in virological sustained responder and non-responder patients

Abstract Background About 130 million people are infected with the hepatitis C virus (HCV) worldwide, but effective treatment options are not yet available. One of the most promising targets for antiviral therapy is nonstructural protein 3 (NS3). To identify possible changes in the structure of NS3 associated with virological sustained response or non-response of patients, a model was constructed for each helicase NS3 protein coding sequence. From this, the goal was to verify the interaction between helicases variants and their ligands. Findings Evidence was found that the NS3 helicase portion of non-responder patients contained substitutions in its ATP and RNA binding sites. K210E substitution can cause an imbalance in the distribution of loads, leading to a decrease in the number of ligations between the essential amino acids required for the hydrolysis of ATP. W501R substitution causes an imbalance in the distribution of loads, leading and forcing the RNA to interact with the amino acid Thr269, but not preventing binding of ribavirin inhibitor. Conclusions Useful information is provided on the genetic profiling of the HCV genotype 3, specifically the coding region of the NS3 protein, improving our understanding of the viral genome and the regions of its protein catalytic site.

Non-classical HLA-E gene variability in Brazilians: a nearly invariable locus surrounded by the most variable genes in the human genome

The non-classical human leukocyte antigen (HLA) class I genes present a very low rate of variation. So far, only 10 HLA-E alleles encoding three proteins have been described, but only two are frequently found in worldwide populations. Because of its historical background, Brazilians are very suitable for population genetic studies. Therefore, 104 bone marrow donors from Brazil were evaluated for HLA-E exons 14. Seven variation sites were found, including two known single nucleotide polymorphisms (SNPs) at positions +424 and +756 and five new SNPs at positions +170 (intron 1), +1294 (intron 3), +1625, +1645 and +1857 (exon 4). Haplotyping analysis did show eight haplotypes, three of them known as E*01:01:01, E*01:03:01 and E*01:03:02:01 and five HLA-E new alleles that carry the new variation sites. The HLA-E*01:01:01 allele was the predominant haplotype (62.50%), followed by E*01:03:02:01 (24.52%). Selective neutrality tests have disclosed an interesting pattern of selective pressures in which balancing selection is probably shaping allele frequency distributions at an SNP at exon 3 (codon 107), sequence diversity at exon 4 and the non-coding regions is facing significant purifying pressure. Even in an admixed population such as the Brazilian one, the HLA-E locus is very conserved, presenting few polymorphic SNPs in the coding region.

MicroRNA miR-146b-5p regulates signal transduction of TGF-beta by repressing SMAD4 in thyroid cancer

MicroRNAs (miRNA) are small non-coding RNAs involved in post-transcriptional gene regulation that have crucial roles in several types of tumors, including papillary thyroid carcinoma (PTC). miR-146b-5p is overexpressed in PTCs and is regarded as a relevant diagnostic marker for this type of cancer. A computational search revealed that miR-146b-5p putatively binds to the 3' untranslated region (UTR) of SMAD4, an important member of the transforming growth factor beta (TGF-beta) signaling pathway. The TGF-beta pathway is a negative regulator of thyroid follicular cell growth, and the mechanism by which thyroid cancer cells evade its inhibitory signal remains unclear. We questioned whether the modulation of the TGF-beta pathway by miR-146b-5p can contribute to thyroid tumorigenesis. Luciferase reporter assay confirmed the direct binding of miR-146b-5p on the SMAD4 3'UTR. Specific inhibition of miR-146b-5p with a locked nucleic acid-modified anti-miR-146b oligonucleotide significantly increased SMAD4 levels in the human papillary carcinoma cell lines, TPC-1 and BCPAP. Moreover, suppression of miR-146b-5p increased the cellular response to the TGF-beta anti-proliferative signal, significantly decreasing the proliferation rate. The overexpression of miR-146b-5p in normal rat follicular PCCL3 cells decreased SMAD4 levels and disrupted TGF-beta signal transduction. MiR-146b-5p overexpression in PCCL3 cells also significantly increased cell proliferation in the absence of thyroid-stimulating hormone and conferred resistance to TGF-beta-mediated cell-cycle arrest. Additionally, the activation of thyroid most common oncogenes RET/PTC3 and BRAF in PCCL3 cells upregulated miR-146b-5p expression. Our results confirm the oncogenic role of miR-146b-5p in thyroid follicular cells and contribute to knowledge regarding the modulation of TGF-beta signal transduction by miRNAs in PTCs. Oncogene (2012) 31, 1910-1922; doi:10.1038/onc.2011.381; published online 29 August 2011

Dengue virus type 3 adaptive changes during epidemics in São Jose de Rio Preto, Brazil, 2006–2007

Global dengue virus spread in tropical and sub-tropical regions has become a major international public health concern. It is evident that DENV genetic diversity plays a significant role in the immunopathology of the disease and that the identification of polymorphisms associated with adaptive responses is important for vaccine development. The investigation of naturally occurring genomic variants may play an important role in the comprehension of different adaptive strategies used by these mutants to evade the human immune system. In order to elucidate this role we sequenced the complete polyprotein-coding region of thirty-three DENV-3 isolates to characterize variants circulating under high endemicity in the city of São José de Rio Preto, Brazil, during the onset of the 2006-07 epidemic. By inferring the evolutionary history on a local-scale and estimating rates of synonymous (dS) and nonsynonimous (dN) substitutions, we have documented at least two different introductions of DENV-3 into the city and detected 10 polymorphic codon sites under significant positive selection (dN/dS > 1) and 8 under significant purifying selection (dN/dS < 1). We found several polymorphic amino acid coding sites in the envelope (15), NS1 (17), NS2A (11), and NS5 (24) genes, which suggests that these genes may be experiencing relatively recent adaptive changes. Furthermore, some polymorphisms correlated with changes in the immunogenicity of several epitopes. Our study highlights the existence of significant and informative DENV variability at the spatio-temporal scale of an urban outbreak.

Molecular evolution of a malaria resistance gene (DARC) in primates

Genes involved in host-pathogen interactions are often strongly affected by positive natural selection. The Duffy antigen, coded by the Duffy antigen receptor for chemokines (DARC) gene, serves as a receptor for Plasmodium vivax in humans and for Plasmodium knowlesi in some nonhuman primates. In the majority of sub-Saharan Africans, a nucleic acid variant in GATA-1 of the gene promoter is responsible for the nonexpression of the Duffy antigen on red blood cells and consequently resistance to invasion by P. vivax. The Duffy antigen also acts as a receptor for chemokines and is expressed in red blood cells and many other tissues of the body. Because of this dual role, we sequenced a 3,000-bp region encompassing the entire DARC gene as well as part of its 5' and 3' flanking regions in a phylogenetic sample of primates and used statistical methods to evaluate the nature of selection pressures acting on the gene during its evolution. We analyzed both coding and regulatory regions of the DARC gene. The regulatory analysis showed accelerated rates of substitution at several sites near known motifs. Our tests of positive selection in the coding region using maximum likelihood by branch sites and maximum likelihood by codon sites did not yield statistically significant evidence for the action of positive selection. However, the maximum likelihood test in which the gene was subdivided into different structural regions showed that the known binding region for P. vivax/P. knowlesi is under very different selective pressures than the remainder of the gene. In fact, most of the gene appears to be under strong purifying selection, but this is not evident in the binding region. We suggest that the binding region is under the influence of two opposing selective pressures, positive selection possibly exerted by the parasite and purifying selection exerted by chemokines.

GH-Releasing Hormone Receptor Gene: A Novel Splice-Disrupting Mutation and Study of Founder Effects

Background: Mutations in GH-releasing hormone receptor gene (GHRHR) are emerging as the most common cause of autosomal recessive isolated GH deficiency (IGHD). Objective: To search for GHRHR mutations in patients with familial or sporadic IGHD and to investigate founder effects in recurring mutations. Methods: The coding region of GHRHR was entirely amplified and sequenced from DNA of 18 patients with IGHD (16 unrelated) with topic posterior pituitary lobe on MRI. Haplotypes containing promoter SNPs and microsatellites flanking GHRHR were analyzed in patients with c.57+1G>A (IVS1+1G>A) mutation of our previously published kindred and also a Brazilian patient and 2 previously reported Japanese sisters with c. 1146G>A (p.E382E) mutation. Results: A novel homozygous intronic GHRHR c.752-1G>A (IVS7-1G>A) mutation, predicting loss of the constitutive splice acceptor site, was identified in two siblings with IGHD. A compound heterozygous c.[57+1G>A];[1146G>A] and a heterozygous c.527C>T (p.A176V) were found in two sporadic cases. Haplotype analysis provided evidence for a founder effect for the c.57+1G>A mutation and independent recurrence for the c.1146G>A mutation. Conclusion: We report a novel splice-disrupting mutation in GHRHR in 2 siblings and provide evidence that all c.57+1G>A (IVS1+1G>A) mutant chromosomes have the same haplotype ancestor, indicating the occurrence of a founder effect in Brazilian patients with IGHD. Copyright (C) 2012 S. Karger AG, Basel

The genetic diversity of Plasmodium malariae and Plasmodium brasilianum from human, simian and mosquito hosts in Brazil

Plasmodium malariae is a protozoan parasite that causes malaria in humans and is genetically indistinguishable from Plasmodium brasilianum, a parasite infecting New World monkeys in Central and South America. P. malariae has a wide and patchy global distribution in tropical and subtropical regions, being found in South America, Asia, and Africa. However, little is known regarding the genetics of these parasites and the similarity between them could be because until now there are only a very few genomic sequences available from simian Plasmodium species. This study presents the first molecular epidemiological data for P. malariae and P. brasilianum from Brazil obtained from different hosts and uses them to explore the genetic diversity in relation to geographical origin and hosts. By using microsatellite genotyping, we discovered that of the 14 human samples obtained from areas of the Atlantic forest, 5 different multilocus genotypes were recorded, while in a sample from an infected mosquito from the same region a different haplotype was found. We also analyzed the longitudinal change of circulating plasmodial genetic profile in two untreated non-symptomatic patients during a 12-months interval. The circulating genotypes in the two samples from the same patient presented nearly identical multilocus haplotypes (differing by a single locus). The more frequent haplotype persisted for almost 3 years in the human population. The allele Pm09-299 described previously as a genetic marker for South American P. malariae was not found in our samples. Of the 3 non-human primate samples from the Amazon Region, 3 different multilocus genotypes were recorded indicating a greater diversity among isolates of P. brasilianum compared to P. malariae and thus, P. malariae might in fact derive from P. brasilianum as has been proposed in recent studies. Taken together, our data show that based on the microsatellite data there is a relatively restricted polymorphism of P. malariae parasites as opposed to other geographic locations. (c) 2012 Elsevier B.V. All rights reserved.